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SUMMARY 


L 

1 The mission requirements for some satellites require that they 

5 spin continuously and at the same time maintain a precise direction 

1 of the spin axis. An analog-computer study has been made of an atti- 

9 tude control system which is suitable for such a satellite. The con- 

trol system provides the necessary attitude control through the use of 
a spinning wheel, which will provide precession torques, commanded by 
an automatic closed-loop servomechanism system. The sensors used in 
the control loop are rate gyroscopes for damping of any wobble motion 
and a sun seeker for attitude control. The results of the study show 
that the controller can eliminate the wobble motion of the satellite 
resulting from a rectangular pulse moment disturbance and then return 
the spin axis to the reference space axis. The motion is damped to 
half amplitude in less than one cycle of the wobble motion. The 
controller can also reduce the motion resulting from a step change 
in product of inertia both by causing the new principal axis to be 
steadily alined with the spin vector and by reducing the cone angle 
generated by the reference body axis. These methods will reduce the 
motion whether the satellite is a disk, sphere, or rod configuration. 


INTRODUCTION 


The mission requirements of some satellites can best be satisfied 
if they spin continuously. An example is a manned space station which 
spins in order to provide the occupants with a simulated gravity field. 

It has been shown that if such a spinning body is disturbed it will 
wobble with an undamped motion. (See refs. 1 to b . ) Therefore, such 
spinning satellites require an attitude control system. For example, if 
the manned space station is equipped with an auxiliary power unit which 
uses a parabolic solar collector, it will be necessary for this collector 
to point continuously at the sun with fairly good accuracy. In order to 
provide control for such continuous pointing, the control system should 
damp the wobble which will result from any outside disturbance and from 
internal movements that result in changes in product of inertia of the 
station with respect to the body axes. It will also be necessary to 
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provide a means of erecting the station to ’die sunline initially and 
each time the station comes out of the earth's shadow. 

Presented herein are results of an analog-computer study of a 
control system which can satisfy these control requirements. The 
control torques are derived from the precession torques generated 
by a spinning wheel. It is assumed that the control of this wheel 
is automatic through the use of closed-loop control, with the command 
information supplied by rate gyroscopes for rate control and a sun 
seeker for control of the position of the satellite axis. 


SYMBOLS 


X,Y,Z body axes 

Xi,Y^,Zj_ axes fixed in inertial space 


X',Y r ,Z' principal axes; identical to body axes when there is no 
product of inertia 

I-y * I v « 1 17 moments of inertia about body X-, Y-, and Z-axis, respectively, 
X Y Z slug-ffcS 

O 

^XY^XZ^YZ Products of inertia, slug-ft 


L 

1 

5 

1 
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T X ,T Y ,T Z body-axis control torques about X**, Y-, and Z-axis, respec- 
tively, ft -lb 


M 


disturbance moment, ft -lb 


p,q,r rates of rotation about X-, Y-, and Z-axis, respectively, 
radians/sec 


co^,o body-axis components of total rotational rate of control- 
wheel angular -momentum vector, radians/sec 


H 

H f 




angular momentum due to spin of control wheel, ft-lb-sec 
constant defined by equations (38] , slug-ft^ 

Euler angles, radians 

constant vector on Xj_-axis (denotes only magnitude when bar 
is removed), used with no dimensions specified 
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ly m,n components of rj along X-, Y-, and Z-axic, respectively, used 

with no dimensions specified 

&Z outer gimbal deflection, radians unless otherwise specified 

5y inner gimbal deflection, radians unless otherwise specified 

Kp control gain, radians /radian/ sec 

Kp 1 constant defined by equations (38) 

L 

1 K2 control gain, radians/unspecified dimension 

5 

1 s Laplacian variable, per sec 

9 

t time, sec 

Q, frequency, radians/sec 

T l/2 time to damp to half amplitude, sec 

P period, sec 

c angle between body axes and principal axes, radians 

Subscripts: 

o initial condition 

X,Y,Z component in X-, Y-, or Z-axis directions 

X*Y f principal axes 

A dot over a quantity indicates differentiation with respect to time, 

DESCRIPTION OF CONTROL SYSTEM 


General Description 

In order to facilitate a proper visualization of the control system 
and its operation, a general description will be given followed by a 
more exact definition with equations. The torque used by this control 
system is the precession torque produced by a spinning wheel. Precession 
torques arise when a spinning wheel is forced to rotate about an axis 
other than its spin axis. 
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A spinning satellite provides a natural situation for the applica- 
tion of a spinning wheel to provide precession torques. In this case the 
spinning of the satellite provides a constantly available angular-rate 
vector which can be combined with the angular -momentum vector of the con- 
trol wheel to provide a continuous torque. The operation of the assumed 


mechanism is as follows. When no torque Is required, the control -wheel 
angular -momentum vector is alined with the spin vector of the satellite. 

When torque is required about a particular body axis of the satellite, 
the control wheel is rotated on an axis parallel to that body axis, 
thereby a component of the control -wheel an$:ular -momentum vector Is 
produced along an axis that is perpendiculai* to both the satellite spin L 

vector and the satellite body axis for which the torque is required. 1 

This situation produces the desired torque, which can be expressed as 5 

a cross product of the satellite spin rate and the control -wheel angular 1 

momentum. This torque is a nearly proportional function of the tilt of 9 


the control wheel, with a constant spin rate of the station and a constant 
angular momentum of the control wheel assumed. 

A sketch of the mechanism studied in this investigation is shown in 
figure 1. In this case a single control wheel mounted in a double gimbal 
is assumed. This type of gimbal mounting w:.ll produce a conflict between 
the operation of control about one axis and that about the other - that 
is, if the inner gimbal is deflected 90 °, rotation of the outer gimbal 
will not change the direction of the control-wheel momentum vector. For 
small deflections of the inner gimbal, the effect is less pronounced. 

It would also be possible to use two control wheels mounted in single, 
mutually perpendicular gimbals and eliminate this conflict. However, 
the single-wheel configuration offers we ight- saving possibilities. 

In this study it is assumed that the control -wheel gimbal angles 
are commanded by signals from rate gyroscopes mounted on the body axes 
of the satellite and by signals from a sun seeker mounted on the spin 
axis of the satellite and rigidly attached to the satellite. The rate- 
gyroscope signals provide damping, and the sun seeker provides the 
necessary signal for alining the satellite with the sunline. 

Equations of Motion 

The general equations of motion are sirrplified by making the following 
assumptions: 

no coupling exists between the force equations and the moment 
equations 


I XZ " X YZ - 0 
% = M z = 0 
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the product of inertia terms IxyP 
Then the rotational equations of motion 


and Ixy 1 ! can be neglected 
for a rigid body are: 


h P + ( x z - + = T x C 1 ) 

I Y^ + {h ~ I z) rp “ I XY q - r + % = % (2) 

I Z r + ( I Y " ^P' 1 “ J XY p2 + I XY q2 = T Z (5) 


The problems studied herein include the effect of introducing a 
rectangular pulse disturbance about the Y-axis and the effects of prod- 
ucts of inertia on the motion of the satellite. Products of inertia 
exist when the X, Y, and Z body axes are not the principal axes of 
the body. For convenience and clarity in the preceding equations only 
an Ixy product of inertia is included, although Ixz and I^y prod- 
ucts could also exist. 

Presumably, a manned space station such as is being considered 
herein will be built so that either the maximum or minimum principal 
axis of inertia will be parallel to the center line of the solar col- 
lector. The movement of an occupant will therefore cause the product 
of inertia with respect to these body axes to change from zero to some 
finite value. In order to investigate this situation, tests were made 
in which a product of inertia I^y was introduced. The movement of 
an occupant would also cause small changes in the body-axis moments of 
inertia, but these effects are neglected herein because they would have 
only a small effect on the results. 

In all the cases considered, it is assumed that initially the 
satellite is spinning about the X-axis, so that the initial spin rate 
is equal to p Q , and that It is the X-axis that should point in the 
reference direction. 


Control Equations 

The precession torque generated by a spinning wheel is equal to 
the product of the moment of inertia of the wheel, its spin velocity, 
and the rate -of -rot at ion components which are perpendicular to the spin 
axis of the wheel. For this investigation the moment of inertia of the 
control wheel and Its spin velocity are assumed to be constant, and only 
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their product, known as the angular momentum cf the control wheel H, 
will be specified. It is assumed that the angular momentum of the 
controller is due entirely to the spin of the control wheel, that is, 
angular momentum due to rotation of the controller on axes perpendicular 
to the wheel spin axis is neglected. As described before, this momentum 
vector of the control wheel will be alined with whichever body axis is 
prescribed, and therefore it is the body-axis components of the control- 
wheel momentum, %, Hy, and H^, that are of concern. The components 
of the precession torque can be assumed equal to the product of these 
angular -momentum components and the orthogona] rates to which they are 
subjected. The equations for the body-axis components of the precession 
torque are 


T X = " H Z U Y 

(V) 

Ty = - HyO^ 

(5) 

T Z = " H y a5 x 

(6) 


The rate terms o^-, and which are the X, Y, and Z components 

of the total rotation of the control -wheel an^^ular -momentum vector, are 
given by the equations 


a)y = p - sin 

(7) 

a>y = q + Sy cos 6^ 

(8) 

Ofc = r + 5 Z 

(9) 


In order to derive the gimbal rates, &y and- first it will be 

necessary to define the gimbal -angle equations. 

The purpose of the control system is to supply damping torques 
and torques which will aline the satellite X-axis with the reference sun- 
line. The torques which are desired are torques about the Y- and Z-axis. 
Since o^r will be a large and steady value as compared with toy and 
<X£, the Y and Z torques can be obtained by calling for momentum 
components and Hy, respectively, which, it can be seen, are 

multiplied by These momentum components are obtained by commanding 


VO H VJ) h tr 1 



7 


a proper orientation of the gimbals carrying the control wheel. The 
proper orientation of the gimbals is a function of rate -gyro scope and 
sun-seeker signals. 

A practical consideration of the output signals from rate gyroscopes 
leads to the conclusion that, for analytical purposes, these signals can 
be considered equal to the body-axis rates q and r. The analytical 
representation of the sun-seeker signals is as follows; A constant 
vector T), coincident with light rays from the sun, is assumed to exist 
on the reference X^-axis. Body-axis components of this constant vector 
can be defined by the equations: 


l = (cos 0 cos 
m = (cos \|r sin 0 sin (j - 
n = (cos \|r sin 0 cos (j) + 


t)T] 

sin \|r cos 0 )tj 
sin \|r sin 0) T 


(10) 

(id 

(12) 


where , 9 , and 0 are Euler angles, taken in that order. A sketch 

showing these components is given in figure 2. The components m and 
n are used to represent the sun-seeker signals. The fact that the 
quantities m and n have the same characteristics that sun- seeker 
signals will have is illustrated by the following two special examples. 

If the satellite is spinning about its X-axis but with the X-axis displaced 
from the reference Xp-axis, then m and n will oscillate about zero with 
a frequency equal to the spin frequency and with a peak amplitude equal to 
the displacement. If the satellite is spinning with the total resultant 
rotation vector on the reference Xp-axis but with the X-.^is displaced from 

the reference line , then m and n will be constant ••> '_th '[m^ + equal 
to the half angle of the cone generated by the X-axis. 

It is not necessary to determine the Euler angles in order to deter- 
mine l , m, and n. These components can be determined from the following 
relationships: 


mr - nq 

( 15 ) 

np - Zr 

( 14 ) 


A = Zq - mp 


( 15 ) 
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l = r| + 



l dt 


(16) 


m = 


m dt 


(17) 


n = J' n dt (l8) 

L 

It is now assumed that the desired damping aid attitude control torques 1 

will be obtained if the gimbal deflections are defined by the control 5 

equations ^ 

5 Z = K x r + K 2 m (19) 


6y = K-^q. - K 2 n (20) 


where Ky and K 2 are the constant control gains. Kp has the dimen- 
sions of radians per radian per sec, whereas K 2 is a nondimensional 
number that, for small deflections, can be thought of as expressing 
radians (of gimbal angle) per radian (of angular displacement of body 
axes from reference space axes) . The mathematical signs given in these 
control equations are for the situation where the control -wheel angular- 
momentum vector is nominally in the same direction as the spin vector of 
the satellite, as shown in figure 1. The ccntrol should be such that the 
gimbal always runs ahead of the satellite, lhat is, the gimbal should tilt 
further than the satellite. 

The body-axis components of the control -wheel angular -momentum 
vector can now be defined as 


Hy- = H cos &y cos *'z 
Hy = H cos by sin 
Hrz — — H sin by 


(21) 

(22) 

(23) 


Also, the rate factors, o^-, o>y, and that are a part of the 

precession-torque terms can now be defined. As was stated previously, 
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these rate factors are made up of the body-axis rates and the gimbal 
rates, which are added vectorially . Since the gimbal angles are func- 
tions of body-axis rates and m and n, the gimbal rates are functions 
of the derivatives of these factors: 


L 

1 

5 

l 
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&z = K^r + (2*0 

&Y = K-ji - Kgn (25) 

Substituting equations (24) and (25) into equations (7), (8), and (9) 
gives 


OJjj- = p - | K 2_4 - k 2^) sin^r + K^m) 

(26) 

ay = q + ^K-]_q - cos^K-jjr + Kgm^ 

(27) 

ar7 = r + Kyf + K-pm 

(28) 


Substituting the expressions for the body-axis components of the preces- 
sion torque (eqs. (4) to (6)) into equations (l) to (3).> and then making 
the substitutions for gimbal deflections and angular -momentum vector 
from equations (19) to (23) and the substitutions for the rate factors 
from equations (2 6) to (28) yields 


I x p + (l^ - Iy)v + Iyypr = H^cos^K^q - KgnJsln^K^r + K^)] + Kfr + Kgi ^ 

+ H^sin^q - K 2 n )j(q. + (Kj4 - K^cos^r + (29) 


Iyq + (i x - Izi)rp - lyyqr + My = - H^Bln^q - Kgnj J |p -(k l 4 - K^njsin^Kyr + KginjJ 

- H^cos^Kyq - KgnjcoB^Kjr + K^mjJ + ICf" + ) (30) 

lyi r + (l Y - Iy) pq - Iyyp 2 + Iyyq 2 = H^coe^q - K 2 n)cos^K 1 r + + (^4 - K^cos^r 

+ Kgm^J - H^cos^Kyq - Kgn) sln^Kyr + K^Jj^p - (lC-,4 - sln^r + Kginjj (31) 
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The original intention was that, for example, a body rate q would 
caul for a Y torque. Generally, this intention is carried out, hut the 
preceding equations show that other torques are also commanded as a result 
of cross -coupling effects. These cross -coupling effects mak e it difficult 
to predict the effect of a given control commard and make it desirable to 
conduct an analytical study of the system. 

The use of simple constants Kq_ and K 2 for the control gains means 

that perfect servomechanism operation is being assumed. Since the highest 
frequencies that will appear in the solutions are equal to the spin fre- 
quency, and since servomechanisms can be made with natural frequencies L 

much higher than this spin frequency, it is felt that this assumption is 1 

adequate for predicting the operation of an actual system. 5 

1 

9 

Linear Equations 

It is possible to write linear equations for a system restricted 
to the principal body axes and the controller used as a damper by making 
some simplifying assumptions. A general solution of the simplified equa- 
tions gives the natural frequency and damping of the system. The neces- 
sary assumptions are: 


p = constant 


XY 


= 0 


k 2 = 0 


cos 


- 1 

si n (K}.q) = K 1<1 


COS 


■M - 1 

sin^K^rj = K^r 


p » (8^4) (Kjr) 
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The equations (30) and 

+ (*x - 

h* + (h - 


( 31 ) then become 

I z )rp = -HK^qp - H(r + I^r^ 

I x ) pq = H(q + K x qj -HK-^rp 


(32) 

(33) 


L 

1 

5 

l 
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The solutions of these equations for an uncontrolled body have been 
treated in other studies, for example, references 1 and 2, but will be 
repeated herein to facilitate comparison with solutions for a controlled 
body. The characteristic equation for the uncontrolled body is 


s 2 + ih. : Stlfcf :* x ) P 2 = 0 w 

3 Tf I Z 

which, for lx > Iy> Iz or < ly^Z defines an undamped oscillation 
with a period given by the expression 


0. 





(35) 


These solutions can be nondimensionalized by differentiating with respect 
to nondimen sional time pt instead of dimensional time t. The results 
will then be expressed in terms of spin revolutions instead of seconds. 

The characteristic equation for the controlled body is 


s 2 + 


2I x HK 1 p + 2H 2 K 1 


hh 




■ jh - - I x)? 2 + + ( 2I X - *Z - + _^ =0 

I Y I Z + 


(36) 


The time to damp to half amplitude is given by the expression 


[ l/2 


_ 0.692(1^ + H%i 2 ) 


2I x HK 1 p 


+ 2H^K, 


(37) 
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The solutions for the controlled body can be nondimensionalized by 
making these additional substitutions: 


H = pH' 


K-i ' i 

k i=-T' 




( 38 ) 


The period and time to half amplitude would then be expressed as revolu- 
tions instead of seconds, and the values obtained for a given configura- 
tion would apply for any initial spin velocity. 


RESULTS 


L 

1 

5 

1 
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Linear Equations 

The period and time to half amplitude obtained by using the linear 
equations are shown in figure 3. The conditions made in these calculations 
are as follows. The rolling velocity p is assumed to be constant at a 
value of 1 radian per second. The curves giver in figure 3 represent the 
results for a variety of axially symmetric configurations. The variation 
in shape is defined in terms of the ratios Iy/Iy or I z /I x . The moment 
of inertia I x is assumed to be constant and the moments of inertia Iy 
and I z are assumed to be equal and to vary sc that the ratios I z /Iy 

and Iy/Iy vary from 0-5 (a flat disk) to 2 (a rod). The control wheel 
is assumed to have an angular momentum equal tc l/20 of the spin momentum 
of the satellite body (^Txp)* This value represents a wheel size and 

wheel spin rate that could easily be carried in a satellite. The gain 
Kp between the body-axis rates q and r and the control -wheel gimbal 
angles & z and 6y is assumed to be 1 radian per radian per second. 

For the uncontrolled body the time to half amplitude is always 
infinite. For the disk configuration the period of the wobble is equal 
to the spin period, 6.28 seconds. As the ratio Iy/Iy is increased 
to 1 (this value represents a sphere) the period increases to infinity. 
Further increase in Iy and I z causes a decrease in period, and this 

result represents the wobble of a rod spinning about its minimum axis 
of inertia. 

The addition of the damper control brings about very short times to 
half amplitude, which vary from 3 seconds to 53 seconds. The longer time 
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L 

1 

5 

l 
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to half amplitude for the la*-ge values of Iy and as compared with 

that for the small values of Iy and Iz is a reflection of the fact 
that under the assumptions made for these calculations the larger values 
of Iy and 1% represent larger bodies, but the momentum of the control 
wheel is held constant. Increasing either the angular momentum of the 
control wheel or the gain of the control will shift the curve to 

smaller values. The period of the characteristic motion with the damper 
is very nearly the same as for the body alone, with the exception that 
the peak asymptote is shifted to a higher value of Iy/Iy. This shift 
reflects the fact that the addition of the control-wheel momentum to the 
spin momentum of the body effectively increases lx a small amount. 


Nonlinear Equations 

Particular solutions for the nonlinear equations of motion were 
obtained by using an analog computer. Three different configurations 
were studied. The results for a nearly spherical satellite 
(iy = 9,500 slug-ft 2 ; Iy = I z = 9,000 slug-ft 2 ), with an initial spin 
rate of 0.6 radian per second, are given in figures 4 to 9- The results 
obtained for a disk configuration (iy = 9,500 slug-ft 2 ; 

Iy = I 7 = ^,750 slug-ft 2 ) spinning on the axis of greatest inertia are 
shown in figure 10. The results for a rod shape (1^ = 9*500 slug-ft } 

Iy = Ig = 14,250 slug-ft 2 ) spinning on the axis of least inertia are 
shown in figure 11. These figures are tracings of analog records. 

In all cases in which a control wheel was included, the angular 
momentum of the control wheel is assumed to be 200 ft-lb-sec. This 
angular momentum could be obtained by using a 32-pound flywheel with a 
7-inch radius turning at 5*000 rpm. The control-system gains were varied 
from 5 to 20 radians per radian per second for the rate signals, and 
0.33 to 3-33 tor the attitude signal in various runs. It should also be 
noted that the ordinate scales used in the figures vary from figure to 
figure. 


Nearly Spherical Configuration 

The response of the uncontrolled nearly spherical body subjected 
to a 10-second rectangular pitching-moment pulse My of 15 foot-pounds 
is shown in figure 4. In these tests in which step moments are used to 
disturb the system, the moments are probably too large to represent any 
particular event that might happen to a satellite, but the large moments 
are used as a severe test of the stability and performance of the system. 
With the use of the linear equations, a wobble period of l 88 seconds is 
predicted for this example, and this prediction is in good agreement with 
the low frequency mode in the particular solution shown. 
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The displacement of the X-axis from the reference X^-axis is equal 

to sin' "^■(\|n 2 m 2 /rj). For convienence in using the analog computer, 

tj was given a value of 100. The magnitude of the maximum values of 
m and n of ±2.6, therefore, indicates angu.ar displacements of ±1.5°. 
The small oscillations superimposed on the l8<$-second wobble oscillation, 
which have a period equal to the spin period of 10.4 seconds, indicate 
that the resultant rotation (spin) vector is displaced from the Xj_-axis 
by a small amount, approximately 0.15°. Such a result would be expected 
as a result of the pitching -moment disturbance. The displacement of the 



The pitch and yaw rates are shown to have max:Lmum values of ±0.0l4 radian 
per second, and therefore the X-axis is displaced from the spin vector by 
an angle of 1.35°* The phase relation of q and r indicates that the 
X-axis displacement from the spin vector alte:~nates from the XY-plane to 
the XZ-plane in quarter periods. A sketch giving a pictorial representa- 
tion of these results is shown in figure 5* 

Next the control system was included (fig. 6(a)), and the system 
was subjected to a pitching -moment pulse of 16 foot-pounds, which this 
time was held on for 36 seconds. The gain on the rate signal in this 
case is 10 radians per radian per second, and the gain on the attitude 
signal is 0.33 radian per radian. The action of the damper control is 
to hold the X-axis close to the spin vector, -thus the body rates r and 
q are kept low. When the disturbance is removed, the X-axis quickly 
alines with the spin vector, and q and r become exactly zero. At 
this time the spin vector is displaced from the X-^-axis by 0.23° as is 
indicated by the ±0.4 oscillations in n and m with the frequency 
of the spin frequency. It should also be pointed o ut that t he varia- 
tions of n and m appear in quadrature, so that 1 'jm 2 + n 2 = ±0.4 also. 

In this case the attitude control is too weak to bring about any reduc- 
tion in this amplitude in the short time of the test. 

An example in which a pitching moment of 120 foot-pounds was applied 
for about 90 seconds is shown in figure 6(b). In this example, the gain 
on the rate signal is 5 radians per radian per second, and on the attitude 
signal is 3*33 radians per radian. When the disturbance is removed the 
body-axis rates and m and n go to zero, and this result indicates that 
both the spin vector and the X-axis are realined with the reference Xi-axis 
It should be noted that in this example gimbaJ. deflections of nearly b0° 
are encountered. Figure 6(c) illustrates more clearly the manner in which 
m and n are brought to zero. Close examination of the analog-computer 
record shows that the short period mode is being slowly damped out. 
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These examples demonstrate that the controller is capable of 
removing the wobble that results from a step disturbance and can realine 
a given body axis with the reference sunline after it has been displaced. 

A disturbance which is likely to occur in a manned space station is 
a movement of the occupant, which will cause a change in product of inertia. 
Presumably, the station would be built so that either the maximum or mini- 
mum principal axis of inertia will coincide with the center line of the 
solar collector. The movement of an occupant would therefore cause the 
product of inertia with respect to the body axes to change from zero to 
some value. Therefore, tests were made in which a product of inertia 
Ixy of 50 slug-ft 2 was introduced. The results for the uncontrolled 
spherical body are shown in figure J. 

With the body-axis moments of inertia of the satellite as given 
(I x = 9,500 slug-ft 2 ; Iy = Iz = 9,000 slug-ft 2 ), the product of inertia 
of 50 slug-ft 2 indicates that the principal axis is displaced 5*73° from 
the X-axis in the XY-plane. The relationship defining this shift is 


I Y - I Y 

= — *■ sin 2c + cos 2c 

With X T and Y’ taken as the principal axes of inertia, Iyiyi is 
zero. Setting ^x’Y' e l ua ^- ^o zero allows the determination of the 

angle c, which is then the angle between the body X- or Y-axis and the 
principal X 1 - or Y T -axis, respectively. 

A spinning body has a preference for spinning on a principal axis 
of inertia. Thus, in the present example, with small factors neglected, 

p 

the yawing torque “lyyP produced by the product of inertia will be 
balanced by the precession torque ( Xy - Iy)pq w ^ en the body yaws 5*73° 

so that the principal axis is under the spin vector. However, when there 
is no damping in the system, an overshoot equal to the initial movement 
will occur. Thus variations in m, which indicate angular displacements 
from 0° to -11. 5°, and variations in q from 0 to -0.12 radian per 
second appear in figure 7* The precession coupling terms bring about 
±5.73° oscillations about the Y-axis, and cause corresponding variations 
in n and r . 

The locus of points on the body which come under the spin vector 
trace a circle around the principal X’-axis. Therefore it can be seen 
that the solution presented herein is in agreement with the development 
known as Poinsot's Construction. (See, for example, ref. 1.) 
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The addition of a damping control tends to bring the principal axis 
into alignment with the spin vector, with the result that the X-axis has 
a steady displacement of 5-73° from the spin vector - that is to say, the 
X-axis describes a 5 * 73 ° cone about the spin Aector. An additional effect 
takes place, however. By a process which wilJ be described, the control- 
wheel angular -momentum vector reaches a steady- state deflection which 
causes a steady-state torque to be applied which opposes the torque pro- 
duced by a product of inertia, and the angle of the generated cone is 
reduced in size. 

A description of this process, which is defined mathematically by 
equations (29) to ( 31 ) , follows. At this point it is convenient to write 
down the significant parts of equations (30) and (31) : 

+ - Iz)rp = -HK-jqp 

+ ^Iy - ly^PQ. “ lyyP = — HK^rp 


The main torque produced by the product of inertia is a torque on 
the Z-axis described by the term IxyP^* Thin torque can be opposed by 

the torque created when an angle 6^ (that in, Kpr) is commanded. 

p 

The relation is - H(K^r)p = the unbalanced torque about the Z-axis. 

In this example, with H = 200 ft -lb -sec and p = 0.6 radian per second, 
a value of 6^ of 8° is required to completely balance the torque produced 
by the product of inertia. 

When the rate gyroscopes alone are used to command the gimbal angles, 
the following process takes place. The product of inertia produces a 
torque about the Z-axis which, in the steady state, causes a pitching 
rate q to be measured. This pitching rate will call for an angle 5y 
which will produce a torque about the Y-axis. This torque will cause a 
yawing velocity to be measured in the steady state. This yawing velocity 
will call for an angle 5^ which will produce a torque about the Z-axis 
which will oppose the product-of-inertia torque and thus reduce the cone 
angle. It is not possible, with the assumed control equations, to reduce 
the cone angle to zero, of course, because the situation is similar to 
that of a weak -position servomechanism acting against a considerable 
load. The possibility exists, however, of including the integral, of the 
attitude error in the control equation and thereby reducing the cone angle 
to zero. 

With a gain of 8 radians per radian per second for the rate signals, 
the steady-state deflection of the X-axis in the XY-plane is 1.6°, and in 
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the XZ-plane is 1.3°> for a resultant cone angle of 2.0k° (fig. 8(a)). 
With a gain of 20 radians per radian per second, the cone angle is 0.7° 
(fig. 8(h)). 

The transient motion brought about by a step increase in Ijry, an, f 

the action of the controller, cause the spin vector to be displaced from 
the reference X-^-axis, and a small oscillation with a period equal to the 
spin period is superimposed on the steady-state values of n and m, as 
is shown in figures 8(a) and 8(b). The addition of the attitude signal 
causes this superimposed oscillation to be slowly eliminated, with this 
elimination indicating that the spin vector is being alined with the 
reference axis. The larger the rate gain, the longer the time required 
to attenuate this oscillation. The attitude signal also causes changes 
in the steady-state gimbal angles, which result in a slightly altered 
cone angle. These results are shown in figure 8(c). 

If the rate -gyroscope input axes were rotated 90° with respect to 
the body axes, so that 


S z == -K x q 

Sy “ K l r 


then the measured effect of a product of inertia would call for gimbal 
deflections such that the control torque would directly oppose the 
disturbance torque. However, there would be no damping in the system. 

A case was tried in which it was assumed that the gyroscopes were advanced 
through a lead angle of 6o°, with a gain of 20 radians per radian per 
second, so that 


6 Z = 20(0. 5r - 0.866q) 


8y = 20(0.5q + 0.866r) 


The results are shown in figure 9- A comparison with the case in which 
the gyroscopes were not rotated shows a reduction in cone angle from 
0*7° to 0.51 5° with the gyroscope rotated. Further comparison shows 
a slight reduction in damping with the gyroscopes rotated. A linear 
analysis predicts such a reduction in damping. Again, the addition of 
the attitude signals reduces the superimposed oscillation to zero. 
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Disk and Rod Configurations 

P 

The motion of a disk, with a product oi inertia I^y of 50 slug-ft , 

is shown in figure 10. In this case the prcduct of inertia of 50 slug-ft 1 2 
indicates that the principal axes are shifted 0.5° from the body axes, and 
therefore the variation in m indicates an oscillation from 0° to -1°, 
and a pitching-rate oscillation from 0 to -0.012 radian per second occurs. 
The wobble frequency is equal to the spin frequency in this case, and the 
combination of wobble motion and spin motion results in a space -oriented 
trace of the X-axis which is a straight line rather than a circle as was 
obtained with the nearly spherical body. This result is indicated by the 
fact that n remains at zero. 

The addition of the controller brings the principal axis into steady 
alinement with the spin vector, and the X-axis moves in a 0-5° cone. The 
controller is not as effective in reducing the cone angle in this case as 
it was in the case of the sphere. The small cone results in small body 
rates q and r, and a very high rate gain would be required to bring 
about the needed to reduce the size of the cone. As can be seen 

in figure 10, the combination of rate signal, and attitude signal results 
in a 52, of 0, and the cone angle is not reduced by any significant 
amount. 

The results obtained with the rod configuration (fig. ll) are very 
similar to those obtained with the disk. In this case, again, the princi- 
pal axis is shifted 0.5° from the body axis by a product of inertia of 
50 slug-ft 2 . With the rod, the wobble period is 32 seconds. The combi- 
nation of the wobble period and the spin period of 10.4 seconds results 
in the peculiar wave shape for m and n. For a given rate gain, the 
time to half amplitude is larger for the rod than for the disk for the 
reasons given in the discussion of the linear equations. 


CONCLUSIONS 


An analytical study of a wide variety of cases of the use of a 
controller, which utilizes precession torques produced by a spinning 
wheel and which is commanded by rate -gyro scope and sun-seeker signals, 
to control a spinning body shows that the controller will perform the 
following functions: 

1. The controller reduces the motion of the body resulting from 

a rectangular pulse moment disturbance by causing the reference body 
axes to be steadily alined with the spin vector, and by realining the 
reference body axes with the space reference line. 
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2. The controller reduces the motion resulting from a step increase 
in product of inertia by bringing the new principal axes into steady 
alinement with the spin vector, and by reducing the cone angle generated 
by the reference body axis. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va. , April 11, 1961 . 
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Figure 1.- Sketch of controller for a spinning body 
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Figure 11.- Motion of a rod with Ixy = 50 slug-ft 





